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Abstract 

The Bohm interpretation of quantum mechanics is apphed to a transmission and 
reflection process in a double potential well. We consider a time dependent peri- 
odic wave function and study the particle trajectories. The average time, eventally 
transmitted particles stay inside the barrier is the average transmission time, which 
can be defined using the causal interpretation. The question remains whether these 
transmission times can be experimentally measured. 
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1 Introduction 

It is well known that (under certain circumstances) the causal interpretation is empir- 
ically equivalent to the orthodox Copenhagen interpretation |l|, |^. However, there are 
a number of physical problems for which this orthodox approach provides no clear-cut 
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answers. One of such problems is the question of a time observable for a tunneling process. 



Tunneling is the quantum mechanical phenomenon that a particle can cross a barrier with 
potential V, even if its energy is strictly less than V. It is a natural question to ask how 
long it takes on average for particles to cross such a barrier. Unfortunately quantum the- 
ory does not provide a suitable time operator, whose expectation value for a given wave 
packet can be compared with experiment. Time enters quantum mechanics as a param- 
eter, not as an operator. Therefore, this question about tunneling time is not an easy one. 



Tunneling processes may be classified in two types: scattering type, where a wave packet 
is initially incident on a barrier, and then partly transmitted; and decay type when, the 
particle is initially in a bound state, surrounded by a barrier, and subsequently leaks out 
of this confinement. Many authors have addressed the duration of tunneling processes, in 
case of scattering processes ^, ^ and in case of a decay process 0. 



In Bohm's causal interpretation of quantum mechanics various concepts of tunneling times 
for scattering processes can be distinguished. The most well-known time is dwell time, i.e. 
the time particles spend inside the barrier. For an ensemble of particles, we can determine 
the average dwell time. This average dwell time {{td)) can be decomposed into an average 
transmission times {{tt)), i.e. the time spent inside the barrier by those particles, which 
eventually cross the barrier, and average reflection time ((tr)), i-e. the time consumed by 
reflected particles 

{U) = \T\'{tt) + \Rf{tr) (1) 

where \Tf and |-R|^ are the transmission and reflection probabilities respectively. This 
relation consists of: 

(f,) = |T|'(ti) + |i?f (V) + |i?"f (V) (2) 

where R' and t^' refer to particles, which penetrate the barrier but re-emerge on the 
same side and R" and t". to particles that do not enter the barrier. This relation assumes 
that particles are either transmitted or reflected, i.e. they do not remain inside the barrier. 



In this paper, we consider a decay type of transmission and reflection by applying the 
causal interpretation to a double potential well. The decomposition of the dwell time 
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into average transmission and reflection time according to relation (|l]) in a decay type 
of tunnelling is not as straightforward as for the scattering cases. Definitions for trans- 
mission and reflection probabilities are not common for this decay type of transmission 
and reflection. However, we shall propose natural definitions for these concepts and use 
these to define average transmission and reflection time. For reflection we will concentrate 
on particles, which penetrate the barrier but re-emerge on the same side. It will again 
appear to be convenient to obtain the average transmission time by using the concept of 
the average arrival time. 

The aim of this paper is firstly, to investigate whether the causal interpretation of quantum 
mechanics provides a straightforward way to define transmission times in double potential 
wells. For simplicity, we will study periodic wave functions. Secondly, the question 
whether it is necessary to adopt the causal interpretation to give meaning to the thus 
obtained transmission times is addressed. The paper is organised as follows. In section 
we apply the causal interpretation method to transmission and reflection in a double 
potential well. In section ^ we define transmission times in terms of the probability density 
of the wave function. The last section is devoted to a discussion about the question 
whether or not the causal interpretation provides a clear way to define transmission times 
in a double potential well and about the necessity to rely on the causal interpretation for 
this definition. 

2 The causal interpretation applied to a double po- 
tential well 

In this section, the causal interpretation is applied to transmission and reflection in a 
double potential well. We discuss the possibihties to define transmission and reflection 
coefficients, average dwell time, average transmission time and reflection time according 
to this method. 

2.1 Description of the double potential well 

To model the decay type of transmission processes we consider a wave packet in a double 
potential well. The potential well is described by a one-dimensional box, defined from —a 



3 



V 



Figure 1: The system under consideration: a double potential well. The total length of 
the box is 2a, the potential at a and —a is infinite. The barrier is situated from— 6 to b 
and has a constant height V. 



to a (see Fig. At these points, the walls are infinitely high. In the middle of the box 
a barrier from —b to b is situated {b < a). 
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if \x\ > a, 
ii b < \x\ < a, 
if < bl < 6 



(3) 



We consider a wave packet which is, initially, concentrated on one side of the barrier only, 
and in the course of time moves to the other side of the wellQ. 



To obtain such a wave packet, we calculate, according to standard procedure (see Ap- 
pendix) the energy eigenvalues and eigenfunctions for this double potential well. A linear 
combination of the lowest even energy eigenfunction {fe{x,t)) and the lowest odd energy 
eigenfunction {fo{x,t)) is the wave function considered in this paper: 

vl/W(x,t) = i=/,(x,t) + i=/,(x,t) (4) 

The choice of these constants assures a high probability of finding the electron at t = 
between —a and —b. 

Eqn. (^) provides a wave packet with (approximately) maximum extinction at one side. 

However, this extinction is not complete and can never be complete. We could enhance 

this extinction by using a linear combination of more than two eigenfunctions. But this 

^ Note that in the causal interpretation a particle described by a stationary wave function does not 
move, and hence no tunneling occurs. 
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might introduce recurrent trajectories, which we want to avoid. 

Our wave packet shows periodic behavior, and transmission from one side to the other 
takes place in half a period. This period time is inversely proportional to E^^^ — 
where and Ef'^ are proportional to the odd and even energy eigenvalues respectively. 
From the wave packet we can calculate trajectories. The particle velocity in the causal 
interpretation is given by H, |^ : 

where the right-hand side denotes the probability current density over the probability 
density function. Analytical expressions of the velocity (i) can be obtained by using the 
expressions given in the Appendix. 

The trajectories are obtained as the solutions to equation (^), subject to the specification 
of the initial condition xq. An ensemble of possible trajectories associated with the same 
wave is generated by varying xq. Hence, to obtain the trajectories x{t) the differential 
equation (Eqn. |]) should be solved after a starting value (xq) is chosen. We solved this 
differential equation numerically [|10]. The trajectories are given in Figs. and |^. In the 
last figure trajectories are weighted by the initial particle density p(x,0) = |^(x,0)p. 

2.2 A classification of the trajectories 

Let us first note some salient aspects of the form of the possible trajectories. Fig. ^ shows 
various trajectories, that are generated by varying the choice of an initial position (xq). 
One clearly sees their periodic behavior, and that the trajectories do not cross. Note that 
at any given time, all trajectories move in the same direction, i.e., to the right during the 
first half period, and towards the left during the second half. The figure shows that at 
very left and at the very right there are trajectories, which remain on the same side of the 
barrier. This occurs for all trajectories with a starting point left of point Si and to the 
right of h. The area indicated with R (returners) contains trajectories, which reach inside 
the barrier but do not pass it, they return to their original area. Their starting positions 
at t = is between si and S2- These trajectories correspond to reflection. The trajectories 
starting in area T (travellers) are trajectories, which pass the barrier and end up, after 
half a period, at the other side of the barrier. At t = the starting point of the trav- 
ellers is between S2 and —h. This behavior corresponds to transmission. The trajectories 



5 



between —b and b are already inside and between b and a already past the barrier at t = 0. 

In Fig. 1^ the probability density |\I'(x,0)p according to the wave function (^) is taken 
into account (the meaning of the symbols a, b, Si, S2 and ti are given in Figs. ^ and 
Along the time axis we marked some special times, tp tm and t„. Time tp is the 
instant at which the trajectory, starting at t = in S2 passes at —b. This trajectory 
marks the bifurcation between the reflection and transmission area. Time tm indicates 
the time the trajectory starting at t = in —6 arrives at b. This is the time needed to 
deliver the trajectories, which at t = already are inside the barrier, to the right hand 
side of the barrier. Although time tn lies outside half a period of time, we marked this 
time because it gives information about reflection times. All trajectories starting at t = 

between si and S2 are at t = ti inside the barrier. Between t = ti and t = tn they leave 

2 2 

the barrier and hence, reflected particles are inside the barrier between time tp and time tn- 



2.3 Definitions of transmission and refiection coefficients 

In view of the above classification, the most straightforward way to define transmission 
and reflection seems to reserve the term transmission for travellers (T) and reflection for 
returners (R). Indeed, only travellers are actively involved in the transmission process. A 
probability for transmission, the transmission coefficient \Tf, can then be defined as: 

\Tf = f ^ dx\^<{x,Q)\^ (6) 
The reflection coefficient can be defined likewise: 

= r rfx|^(x,0)|^ (7) 



Note that this reflection coefficient is in fact related to 



in relation fEI). The trans- 



mission and reflection coefficients (H and ^ do not add up to unity: \T\^ + \R\^ 7^ 1. Of 
course, this is due to the fact that we did not include R and that there is a finite prob- 
ability that the trajectories are already inside or past the barrier. These last possibilities 
are usually excluded in discussions of tunneling in one- dimensional scattering processes 

m. 
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Figure 2: Bohm trajectories to show the different starting point Xq possibiUties at t = 0. 
Left of S2 and right of b trajectories do not leave their own area. R=reflection area, 
starting positions between si and S2- T=transmission area starting positions between S2 
and —b. 
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Figure 3: Bohm trajectories for a one dimensional double potential well. We represent 
the probability density function t = by N points. These points are spaced with equal 
probability density intervals. J^^ dx \^{x, 0)| = -j^- Xn is the position of Xq dependent 
on A^. Ntot is the total number of points. The value of — | is arbitrarily chosen. Any other 
value (between and 1) gives an equally valid description, tp is the time the trajectory, 
starting at t = in S2, passes at —h. tm = ti — tp and tn = ti + tm- In this example 
N = 15. 
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2.4 Definitions of dwell time, transmission and refiection times 

Average dwell time is the average time that particles spend inside the barrier region: 

{td) = dt dx\m{x,t)Y (8) 

Jo J~b 

The probability to encounter a particle inside the barrier {j\ dx |\l'(x, t)^) is for the wave 
function given in Eqn. (|) independent of time and hence, the average dwell time inside 
the barrier in half a period is the probability to encounter a particle inside the barrier 
multiplied by half a period: 

{td) = U t dx\^{x,t)\^ (9) 

2 J —b 



To define average transmission and reflection times, we have to go back to the trajectories 
because trajectories provide information about the position {x) at each instant of time 
(t), which, under the assumption that each trajectory passes x only once, can be inverted 
to give the function t{x). It is convenient to express the transmission and refiection 
times in terms of arrival time distributions. In particular, let the instant at which a 
trajectory, starting in xq at t = 0, arrives at xi, be denoted as t(xo; xi). Averaging over 
the probability density |\I^(xo,0)p that a particle starts at Xq, we obtain the arrival time 
distributions (n(t)): 

_ / rfxo|^(zo,0)p^(t-t(3;o;xi)) 

jdxo\^{x,M' ^ ' 

where the integration limits should be chosen in such a way as to fulfill the assumption 
that each trajectory that passes Xi does so only once in half a period. 
Hence, the arrival time distribution of the to-be-transmitted particles at the exit of the 
barrier (6) is: 

n ^ Is-'dxo\^ixo,0)m-tixo;b)) 
J-' dxo\^{xo,OW 

The denominator is the transmission coefficient (see formula^). The starting points Xq 
are within the transmission area, between —6 to S2. 

Similarly, the arrival time distribution of the to-be-transmitted particles at the entrance 
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of the barrier (—6) is: 

n dxo\<il{xo,0)mt-t{xo;-b)) 

As long as the trajectories cross a particular point (xi) only once, Leavens ^ showed 
that: 

dxo\^{xo,Q)\^5{t{xi) -t{xi\x)) = j{xi,t{xi)) (13) 

which gives us the arrival time distributions in term of probability current densities. 
Hence, the average arrival time at point xi {{ta{xi))) under the same conditions, is: 

In our case, the probability current density in the double potential well is unidirectional 
for half a period and hence, we can determine the average arrival time at the entrance 
and the exit of the barrier. Taking the time boundaries from Fig. ^ in consideration the 
average arrival time at the entrance of the barrier is: 

= (15) 

lo dtj{-b, t) 

and at the exit: 

im) = (16) 
u;.dtnb,t) 

and hence, the average transmission time, the average arrival time at the exit of the 
barrier minus the average arrival time at the entrance of the barrier, reads: 

«,) = (aw-(M-w = 4^«-^;^^^ (IT) 

To determine reflection times, we have to extend the observation time to t„ (see Fig. 
To-be reflected particles enter the barrier between time tp and ti, they pass the barrier 
again on their way back between time ti and t„. Hence, the average reflection time is: 

Jl-dttj{-b,t) dttj{-b,t) , , 

I,,dtj{-b,t) J^^^dtj{-b,t) 
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In this definition we use time t„, which does not fall inside the range defined as half a 
period (from t = to ti). This might seem in contradiction to Eqn (j^), the definition of 
average dwell time, which uses ti as its upper bound limit. However, the time taken by 
the reflecting particles on their way back to their original place {ti to t„) is equal to the 
time taken by particles already inside the barrier at t = 0, to leave the barrier. The last 
mentioned time was not accounted for as belonging to average transmission time. Hence, 
the integration limits of the addition of average transmission and reflection times and the 
total average dwell times are consistent. 

3 Transmission time in terms of probability density 
of the wave function 

In the previous section, the definitions of the average transmission and reflection times, 
using the average arrival time distribution were based on trajectories from the causal in- 
terpretation of quantum mechanics. However, the final expressions (Eqns. and |TB|) do 



no longer depend on the trajectories but on probability current densities. This suggests 
the possibility to define average transmission (and reflection times) without explicit use 
of the trajectories. In this section, we will show, that average transmission times can, 
indeed be obtained without explicit reference to trajectories. Instead we only rely on the 
non-crossing property of the causal interpretation and the unidirectionality between t = 
and t = ti of the current density flow and the periodicity of the wave functions. Aver- 
age reflection times will not be considered, because we are interested in the transmission 
process. The question whether, the fact that the trajectories are not needed explicitly in 
order to determine the transmission times implies that the causal interpretation is super- 
fluous for this purpose is left for the next section. 



In order to show how average transmission times for a periodic wave function can be 
determined from the probability density, we refer to Figs. ^ and ^ 0. Here we have par- 
titioned the interior of the double well in 5 areas, as shown in the figures. In Fig. ^ 
the probability densities at t = and in Fig. |^ at t = ti are given. One can see that 

^Actually, we chose the constants (see Eqn. (^) in such a way that the probabihty to encounter parti- 
cles at i = at the right-hand side of the barrier is large enough to be visible in graphical presentations. 
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Figure 4: Probability density at t = 0. The transmission area (T) and reflection area 
{R, between si and S2) are indicated. Their destinations are given in Fig. ^. —a and a 
are box boundaries and —b and b the barrier boundaries. 
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-a -b b a 

box 

Figure 5: Probability density at t = ti. Compare to Fig. ^: T is past the barrier and R 
inside the barrier. 

the probabihty density goes from left to right during half a period. The probability den- 
sities in Figs. § and 1^ are mirror images. We will exploit this symmetry in our calculations. 

The 5 areas indicated in Figs. ^ and § have the following meaning: 

• The probability in area at the utmost left-hand side, between —a and Si, remains in 
its own domain, si can be found by the condition that the probability between —a 
and Si should be equal to the probability between —a and —b in Fig. ^. Because 
of symmetry the probabihty density between —a and Si is equal to the probability 
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density between b and a at t = (Fig. Hence, 

r dx\^{x,0)f = r dx\^{x,0)f (19) 

J —a Jb 

• The area which is indicated with an R in Fig. ^ is the reflection area (between si 
and 52)- In Fig. ^ the probabihty within that area has moved inside the barrier. 
The probabihty density inside the barrier is constant at all times and hence is the 
same as the probability inside the barrier at t = 0: 

r dx\<iJ{x,0)f = f da;|^(x,0)|^ (20) 

J s\ J —b 



• We indicated a T in Fig. |^ for the transmission area (between S2 and —h). The 
probability density from this area arrives at t = ti, Fig. |, at the opposite side 
of the barrier. Because the probability density of all other areas are known, this 
probability density can be determined by: 

j ^ dx\^{x,Q)\^ = 1-2 r dx\^{x,Q)\^ (21) 

J S2 J —b 

To determine the average transmission time, we need the average arrival time of the to-be 
transmitted particles at the entrance of the barrier and at the exit of the barrier. We use 

the assumption that the flux is unidirectional between t = and t = ti. 

2 

• At the entrance of the barrier the transmission flux starts at t = and should be 
stopped at the time, when the left-hand side of the barrier is emptied of all travelling 
probability. The probability left behind at the left-hand side of the barrier at t = tp 
is equal to the probability between —a and S2 at t = 0. Hence, this time, tp, can be 
found implicitly by: 

/ dx\'^{x,tp)f = r dx\<il{x,0)f (22) 

J —a J —a 

• At the exit of the barrier, we must wait until all the probability initially inside the 
barrier has passed before the travelling part arrives. This happens between time 
t = and t = tm- After half a period the travelling part is inside the right-hand 
side well. Between time tm to ti the transmission part passes at the exit of the 
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barrier, starts when the right hand side of the barrier contains the probabihty 
of the right-hand side and the probabihty under the barrier: 

r dx\^{x,tm)f = r dx\'^{x,0)f (23) 

The right-hand side terms of Eqns. (|22D and (|23D are equaL The half period 
evolution (ti) is ^^^^^ (see the Appendix for symbols) [|. 

Now, Eqn. ([TTD can be filled in, which gives us the average transmission time of the 
double potential well, without using the trajectories. 



4 Discussion and conclusions 

In this section, we discuss the question whether or not the causal interpretation pro- 
vides an unambiguous way to define average transmission times in a double potential 
well. Secondly, we discuss the necessity to adopt the causal interpretation to define this 
transmission time. We also discuss the experimental accessibility of transmission times. 
In the causal interpretation of quantum mechanics, the position of an individual particle 
travelling along a particular trajectory is determined at each instant of time. The tra- 
jectories show the possible behavior of particles inside the barrier and hence, the causal 
interpretation creates the possibility (at least numerically) to discriminate between trajec- 
tories inside the barrier, whose fate is transmission and whose fate is reflection. Although 
definitions for transmission and reflection coefficients are not common practice in a dou- 
ble potential well, the definitions, given in Eqns. (|^) and (|^) provide a useful tool to 
differentiate between these two possibilities. The wave functions were made out of two 
eigenfunctions and hence the trajectories show that the current density flow, between 
t = and if: = ti , is unidirectional. Hence, a straightforward way to deflne transmission is 
provided. We took the lowest energy level eigenfunctions but any pair of eigenfunctions 
would give an equally straightforward way to deflne transmission. However, the addition 
•^Using the probability current density the times tp and tm can also be found by: 

Jo Jt^ 
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of more eigenfunctions to create wavefunctions would cause recurrent trajectories. In 
case of recurrent trajectories the labeling of transmitted and reflected trajectories is not 
straightforward. 

Similar to the more well- known case of scattering processes, the transmission times can 
be conveniently expressed in terms of the arrival time distributions. Average arrival times 
at the entrance and exit of the barrier are used to express the transmission time. Hence, 
the average time spent inside the barrier by eventually transmitted particles can (for the 
above mentioned wave functions) unambiguously be defined and hence the answer to the 
question whether the causal interpretation can give an unambiguous definition of average 
transmission time in a double potential well is "yes". 

Next, let us discuss the question whether or not the causal interpretation of quantum me- 
chanics is needed to define transmission time. Transmission times for a double potential 
well are defined in the causal interpretation of quantum mechanics. However, the causal 
interpretation also implies a different world view than the standard, orthodox interpreta- 
tion. The trajectories describe the way particles move. The initial position of the particle, 
although unknown to us, fixes its future path completely and hence, in contrast to the 
orthodox interpretation, the causal interpretation theory is a deterministic theory. 

In section ^ we showed how the definition of transmission times could be obtained, without 
trajectories, from the probability density of the wave function. This may suggest that the 
causal interpretation is superfluous for the determination of transmission times. However, 
transmission times are obtained under the assumption that the flow of probability density 
follows the non-crossing property of the trajectories of causal interpretation. Outside the 
causal interpretation, the justification for this assumption is not clear. Hence, our discus- 
sion of the definition of average transmission time in a double potential is dependent on 
this aspect of the causal interpretation of quantum mechanics for its justification. 

Finally, we pose the question of the experimental accessibility of average transmission 
times. For a double potential well model, a definition for the average transmission time is 
offered in this paper but the question whether an experimental set-up to measure average 
transmission times can be devised is open and hence, the question remains whether the 
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average transmission time for a double potential well can be verified by experiments. The 
usefulness of the definition of average transmission times would be greatly enhanced if 
average transmission times can be experimentally measured. 
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Appendix 

We take atomic units, i.e. we put m = 1 and h=l. Taking the boundary conditions in 
consideration gives us even (symmetric) and odd (asymmetric) solutions. For the even 
solutions we find: 

{sin(A;e(a; + a)) ii — a < x < —b, 

sin{K{b-a)r^^ if -b<x<b, (25) 

— sin(A;e(a; — a)) if b<x<a 

Similarly, for the odd solutions, one obtains: 

sin(fco(x + a)) if — a < x < —6, 

foix,t) = NJe-'^"' { sm{K{b - a))l^^ if - b < x < b, (26) 

sm{ko{x — a)) if b < x < a 

where a, b and V are explained in Fig |l], N^^o are (complex)normalisation factors, k^ o = 



2Ee,o, 0:e,o = \/2(^ — Ee^o), and Ee^o < V. Further, for even. Eg is the solution of the 
equation: 

arctan(^^^= coth(6y2(l/ - E^)) = rnr - {a - b)^2Ee (27) 
\/V — Ef. 

and for odd, Eo is determined by 



arctan(— =^=^tanh(6^2(y - Eo)) = rnr - {a - b)^2Eo. {2i 
\JV — Eo 
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Figure captions 
Fig. 1 

The double potential well. The total length of the box is 2a, the potential at a and —a is 
infinite. The barrier is situated from— 6 to h and has a constant height V. 
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Fig. 2 



Bohm trajectories showing the different types of behaviour, depending on starting point 
Xq at t = 0. to the left of S2, and right of b, trajectories do not leave their own area. 
R=reflection area, starting positions between si and S2- T=transmission area starting 
positions between S2 and —b. 

Fig. 3 

Bohm trajectories for a one dimensional double potential well. The probability density 
function is represented at t = by N points. These points are spaced with equal proba- 
bility density intervals. J^^ dx |\E'(x,0)| = -j^- xn is the position of xq dependent on 
A^. Ntot is the total number of points. The value of — | is arbitrarily chosen. Any other 
value (between and 1) gives an equally valid description, tp is the time the trajectory, 

starting at t = in S2, passes at —b. tm = ti — t„ and tn = ti + tm- In this example 

2 ^ 2 

A^= 15. 
Fig. 4 

The probability density at t = 0. The transmission area (T) and reflection area (i?, 
between si and S2) are indicated. Their destinations are given in Fig. ^ —a and a are 
box boundaries and —b and b the barrier boundaries. 

Fig. 5 

The probability density aX t = t\. In comparison with Fig. ^: T is now past the barrier 
and R inside the barrier. 



19 



